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We show that, because of the inevitable twist and lattice mismatch in heterobilayers of transition 
metal dichalcogenides, interlayer excitons have six-fold degenerate light cones anomalously located 
at finite velocities on the parabolic energy dispersion. The photon emissions at each light cone 
are elliptically polarized, with major axis locked to the direction of exciton velocity, and helicity 
specified by the valley indices of the electron and the hole. These finite-velocity light cones allow 
unprecedented possibilities to optically inject valley polarization and valley current, and the obser¬ 
vation of both direct and inverse valley Hall effects, by exciting interlayer excitons. Our findings 
suggest potential excitonic circuits with valley functionalities, and unique opportunities to study 
exciton dynamics and condensation phenomena in semiconducting 2D heterostructures. 

PACS numbers: 74.78.Fk, 78.67.Pt, 71.35.-y, 72.25.Fe 


Monolayers of group-VIB transition metal dichalco¬ 
genides (TMDs) have recently emerged as a new class of 
direct-gap semiconductors in the two-dimensional ( 2 D) 
limit HHS]. These hexagonal 2 D crystals have exotic 
properties associated with the valley degeneracy of the 
band edges, including the valley Hall effect mn\, the 
valley magnetic moment IHHII], and the valley optical 
selection rules [i [Mm, leading to rich possibilities 
for valley-based device applications. The visible range 
bandgap further makes these 2 D semiconductors ideal 
platforms for optoelectronics. Due to the strong Coulomb 
interaction, the optical response is dominated by exci¬ 
tons, the hydrogen-like bound state of an electron-hole 
pair. The demonstrated electrostatic tunability and op¬ 
tical controllability of valley configurations of excitons in 
monolayer TMDs have implied new optoelectronic device 
concepts not possible in other material systems [IMZ]. 

Stacking different TMDs monolayers to form van der 
Waals heterostructures opens up a new realm to extend 
their already extraordinary properties [ 18 ]. M0X2/WX2 
(X=Se, S) heterobilayers have been realized [TM 25 ] . 
which feature a type-II band alignment with the con¬ 
duction (valence) band edges located in M0X2 (WX2) 
layer. Exciton then has the lowest energy in an in¬ 
terlayer configuration (i.e. electron and hole in differ¬ 
ent layers), from which luminescence is observed [TM 2 T] . 
Due to the spatially indirect nature, interlayer excitons 
in MoSe2/WSe2 heterobilayers have shown long lifetime 
exceeding nanosecond, repulsive interaction, and electro¬ 
statically tunable resonance m, all of which are highly 
desirable for the realization of excitonic circuits and con¬ 
densation [ 26 f[ 29 ] . An unprecedented aspect of this inter¬ 
layer exciton system is that the heterobilayers in general 
have incommensurate structures due to lattice mismatch 


and twist in the stacking which, together with the valley 
physics inherited from the monolayers, bring in radically 
new properties. 

Here we discover anomalous light coupling properties 
of interlayer excitons in twisted M0X2 /WX2 heterobilay¬ 
ers. We find these excitons have unique six-fold degener¬ 
ate light cones, located at finite center-of-mass velocities 
on the parabolic energy dispersion. In these light cones, 
an interlayer exciton can directly interconvert with an el¬ 
liptically polarized photon (without phonon or impurity 
assistance), with helicity depending on the valley indices 
of the electron and the hole. These light cones allow reso¬ 
nant optical injection of excitonic valley polarization and 
valley current with versatile controllability, as well as the 
observation of valley Hall and inverse valley Hall effects of 
interlayer excitons. Contrary to known exciton systems, 
here the exciton lifetime is long at the energy minimum 
while optical injection and probe are allowed at finite 
velocities, suggesting unique opportunities to study exci¬ 
ton dynamics and condensation phenomena. For lattice 
matching heterobilayers of AA-like ( 0 ° twist) or AB-like 
( 60 ° twist) stacking, all light cones merge into a single 
one, and the quantum interference leads to dependence 
of the brightness and the polarization selection rules on 
the translation between the layers. 

The optically active interlayer excitons are the ones 
with electrons (holes) from the ±K' (±K) valleys in 
M0X2 (WX2), where layer-hybridization is substantially 
quenched by the large band offsets between the layers 
131 [3QII35] . The term light cones here refer to the exciton 
phase space regions where interconversion with photon is 
directly allowed (without phonon or impurity assistance). 
Concerning the light coupling, interlayer exciton in heter¬ 
obilayers is distinct from existing systems by the fact that 
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FIG. 1: (a) With a small twist and/or lattice mismatch, the 
Brillouin zone corners of the M 0 X 2 layer (r'K', open dots) 
shift away from those of the WX 2 layer (tK, solid dots). An 
electron and a hole can recombine if their momentum sum 
k + k' equals the displacement vector between a tK and a 
T K' corners (green arrows), (b) The interlayer Coulomb in¬ 
teraction (wavy line) between the electron and hole conserves 
their momentum sum [40]. An interlayer exciton formed 
with kinematic momentum Q = ±Qo, iCsQo, or ±(7 |Qo 
can directly recombine to emit photon, (c) These light cones 
at finite velocities are illustrated by the bright spots on the 
energy dispersion, while the exciton is optically dark at the 
energy minimum at Q = 0. The dashed arrows illustrate the 
relaxation path ways. 



FIG. 2: (a) With the elliptically polarized selection rules, a 
linearly polarized light inject exciton currents with different 
rates and directions at the six light cones. The net effect 
is a pure valley current (green and pink colors denoting the 
valleys), (b) Valley Hall effect in heterobilayers. With the 
light cones at finite velocities, the spatial pattern of emission 
polarization also depends on the exciton distribution in Q- 
space. (c) Exciton valley current can induce number current 
through the inverse valley Hall effect, observable from the 
spatial map of the luminescence. 


a small twist and/or lattice mismatch shifts the electron 
valleys away from the hole valleys (Fig. [^a)). The ex¬ 
citon at zero velocity is therefore indirect in momentum 
space. The small mismatch between the electron and hole 
valleys can nevertheless be compensated by the exciton’s 
kinematic momentum Q (Fig. [^b)), to allow radiative 
recombination with a direct transition. This happens at 
Q = ±Qo, TCsQo, or TCiQo, where Qo = rK —r'K' is 
the displacement from the M 0 X 2 Brillouin zone (BZ) cor¬ 
ner r'K' to a nearest WX 2 BZ corner rK, with r' = r 
(r' = —r) if the twist angle 0 is near 0° (60°). Thus, 
in the phase space parametrised by Q, the light cones 
are all located at finite velocities with the six-fold de¬ 
generacy (Fig. [^c)). With low symmetry of the exciton 
states, photon emission at each light cone can have a gen¬ 
eral form of being elliptically polarized (circular or linear 
polarization being the two limits when the ellipticity e 
approaches 0 or 1). The exciton states at Qo, CsQo, and 
CfQo have the same electron-hole valley configuration 
(r', r) and are related by the Cs rotation, so their ellipti¬ 
cal polarizations are also Cs rotations of each other. The 
rest three light cones are their time-reversal counterparts 
with the valley configuration (—r', —r) (Fig. ia)). 

Optical injection of valley current .—With the helicity 
determined by the valley index and the major axis of the 
elliptical polarization locked to the exciton moving direc¬ 
tion, the finite-velocity light cones make possible valley 
selective injection of valley polarization as well as valley 
current (Fig. |^a)). Under excitation by a linearly po¬ 
larized light, the absorption rates are all different at the 
six light cones associated with different velocities, giv¬ 
ing rise to valley current (Fig. |^a)). Using 0 to denote 
the angle between the major axis of polarization of the 
light cone at Qo = Qoy and that of a linearly polar¬ 
ized pumping light, a pure valley current is injected with 


i*ate f in the direction x sin 20 + y cos 20, 

a being the population injection rate. 

Optical valley Hall effect .—This exciton system is an 
ideal platform to study the Berry phase effect in the 
Bloch band. TMD monolayers have valley dependent 
Berry curvatures in their conduction and valence band 
edges, which give rise to the valley Hall effect of the car¬ 
riers laiT]. The exciton can inherit the Berry curvature 
of its electron and hole constituents [31113 , giving rise 
to an exciton valley Hall effect that is observable from 
the spatial pattern of emission polarization (Fig. Hb)) 
j36[ [38] . The elliptically polarized emission on the two 
edges with opposite helicity is a phenomenon unique to 
this heterobilayer interlayer exciton system. The optical 
injection of excitonic valley current through the finite- 
velocity light cones further makes possible the observa¬ 
tion of inverse valley Hall effect, where the valley current 
induces a number current in perpendicular direction due 
to the skewed motion by the Berry curvature (Fig. He)). 

We present now the theoretical details that backup the 
above general statements on the anomalous light cones. 
The heterobilayers can be characterized by the twist an¬ 
gle 0° < ^ < 60° and the in-plane translation tq between 
the layers (Fig. |^a)). The origin of the in-plane coordi¬ 
nates is chosen at a Mo atom, and the coordinate of its 
nearest W atom defines tq m- We first give the descrip¬ 
tion of interlayer valley excitons in twisted heterobilayers 
at vanishing interlayer coupling, and then analyze their 
light-interaction effects due to the residue interlayer cou¬ 
pling. Using \er',k'} (|^r,k)) to denote the electron (hole) 
Bloch state at momentum k' (k) measured from r'K' 
(rK) in M 0 X 2 (WX 2 ) layer in the absence of interlayer 
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FIG. 3: (a) A M 0 X 2 /WX 2 heterobilayer can be specified by the twist angle 0 and the interlayer translation ro. (b) The 
bright spots from large to small denote respectively the main, 1st Umklapp, and 2nd Umklapp light cones in the Q-space, for 
WSe 2 /MoSe 2 heterobilayers with 0 near 0° or 60°. (c) Solid (open) dots denote tK + G (r'K' + G') points of WSe 2 (MoSe 2 ) 
layer in the extended zone scheme for 0 = 5 °. The locations of main (nth Umklapp) light cones are given by the displacements 
between nearest solid and open dots on the Co (Cn) circles, (d) Low energy light cones at ^ = 0°, 5°, 22.8°, 37.2°, 55°, 60°, 
respectively. Evolutions of the main (2nd Umklapp) light cones with 0 are shown by the thick (thin) red lines, (e), (f) Lattice 
mismatched MoS 2 /WSe 2 heterobilayer. All light cones are always at finite Q. (g) Recombination pathways of interlayer 

exciton X _ ^ at 0 near 60°. Solid (dashed) arrows denote dipole transition (interlayer hopping). The inset shows the elliptical 

polarizations of emissions at the main light cones, (h) X++ at 0 near 0°. At the finite Q light cones, the polarization and 
strength of the optical dipole are independent of ro. 


hybridization, the exciton wavefunction writes [311101 


AQ 


with the eigenenergy Ei{Q) = Mq = 

rrie + rrih is the exciton mass, Ag the band-gap, Ei, the 
binding energy, and ^/(AQ) describes the electron-hole 
relative motion with AQ = Besides the elec¬ 

tron and hole valley indices r' and r, an interlayer ex¬ 
citon state is labeled by Q = k + k', the kinematical 
momentum associated with the center-of-mass group ve- 
locity (xg’^K^i-e + [40], which is 

a good quantum number even in absence of translational 
invariance at general twist angles. 

The Fourier expansion of the electron (hole) 
Bloch function consists of the plane-wave components 
g2(k +r'K +G )-r j^g2(k-rK-G)-r^^ whcrC G' (G) dcUOtCS 

all possible reciprocal lattice vectors in the M 0 X 2 (WX 2 ) 
layer. Neglecting the small photon wave vector, an 
electron-hole pair |eT-',k')l^r,k) can recombine only when 


k' + r'K' + G' = —k + rK + G. So the exciton light 
cones are at the Q-space points: 

Q = rK-r'K', (2) 

where t'k' = r'K' + G' (tk = tK + G). 

In Fig. [^c) and (e), tk and t'k' are shown in the 
extended zone scheme for a MoSe 2 /WSe 2 heterobilayer 
at = 5° and a MoS 2 /WSe 2 heterobilayer at = 0°. 
When 0 changes, t'k' (tk) move on concentric circles de¬ 
noted as C'^ {Cn)‘ The nearest (tk, t'k') pairs on Co 
and Cq correspond to main light cones. And the near¬ 
est (tk, t'k') pairs on Cn and C^ define the nth Umklapp 
light cones^ where light coupling is an Umklapp type pro¬ 
cess (Fig. [^b)). Fig. I^d) and (f) shows the locations 
of these light cones as a function of twist angle in a 
range Q < O.liF that corresponds to a kinetic energy 
< 60 meV. The six main light cones are located at 
Qo = (i^ — K' cos 50)^ — K' sin 50y and its ^ rotations, 
where 50 is the deviation of 0 from 0° or 60°. Three of 
them correspond to the same valley configurations (red 
lines in Fig. [^d) and (f)), and their time-reversals give 




















































4 


the remaining three. Ist-Umklapp light cones also have 
six-fold degeneracy, located at 2Qo. The main and 1st- 
Umklapp light cones are the lowest energy ones at 0 near 
0° or 60°. There are twelve 2nd-Umklapp light cones 
(Fig. [^b)), and six of them become the lowest energy 
ones for lattice matching heterobilayers when 0 is near 
the commensurate angle 21.8° or 38.2° (Fig. |^d)). 

At these light cones, the exciton acquires a finite op¬ 
tical transition dipole from the interlayer coupling. The 
interlayer hoppings (Ht) of electron and hole will correct 
the interlayer exciton wavefunction in Eq. 0, hybridiz¬ 
ing q) with small components of the intralayer val¬ 
ley excitons m- The optical transition matrix element 
of the interlayer exciton is then: 


V,,r,Q =(0|D|x(?tQ) + (0|D|X.,) 


-E'/(Q) — Em 


(0|D|X.) 


(X.|gT|4v.Q) 

Ei{Q) - Ew 


( 3 ) 


Here D is the electric dipole operator and |0) is the vac¬ 
uum state. The first term on the right-hand-side is from 
the transition dipole between le^-z^k') and which 

is small because of the spatial separation. The second 
and third terms describe respectively the light coupling 
mediated by bright intralayer exciton in M 0 X 2 layer 
through hole interlayer hopping, and X^- in WX 2 through 
electron hopping. Em and Ew are the energies of the 
intralayer excitons, a few hundred meV above Ej as de¬ 
termined from optical and scanning tunneling measure¬ 
ments [T9H2TJ[30]. can be extracted from ab initio 

calculations at several commensurate stackings [40] , from 
which we can determine the dipole strength, and the el- 
lipticity and major axes of the emission polarization of 
the six-fold light cones at a general twist angle. 

In the two-center approximation [48[ |49| : 


(0|D|x(?tQ) oc 

K, k' 

K, k ' 


where t^^, are the Fourier components of hopping 

integral (transition dipole) between a W and a Mo d- 
orbitals with magnetic quantum numbers m and m' re¬ 
spectively, which decay fast with the increase of |k| [40] . 

Fig. gh) schematically shows interlayer excitons X^-'r 
with valley configuration r' = r. These are the ones in 
the main and 1st Umklapp light cones at 0 near 0° and 
in 2nd Umklapp light cones at 0 near 38.2°. Fig. ig) 
shows X^M with r' = —r, which are the ones in the 
main and 1st Umklapp light cones at 0 near 60° and in 
2nd Umklapp light cones at 0 near 21.8°. The emission 



FIG. 4: Dependence of the transition dipole strength and 
polarization as a function of layer translation ro in lattice 
matching MoSe 2 /WSe 2 heterobilayers, (a) X++ at AA-like 
stacking {0 = 0°). The three panels show respectively the 
dipole strength, the degree of polarization in circular basis, 
and the ellipticity and major axis (representd by short lines) 
of elliptical polarization, from ah initio calculations [54]. (b) 
Same plots for X _^ at AB-like stacking {0 — 60°). 


polarizations of the main light cones are illustrated in 
the insets. For MoSe 2 /WSe 2 heterobilayers, we find the 
emission polarization of X++ at 0 near 0° has an elliptic¬ 
ity e = 0.35, which correspond to a ratio of 5 between the 
e+ = and e_ = component of The 

e+ component is contributed by all three terms in Eq. 
0, while the much smaller e_ component is from the 

(0|D|x|?| q) term only. X _is simply the time reversal 

of X++. At 0 near 60°, the valley configurations _ 

and X _^ have the ellipticity e = 0.33. 

Erom ab initio calculations m, X±± at 0 near 0° has 
a transition dipole strength of ^ 0.052D at the main light 
cones, and ^ O.OOII) at both the 1st and 2nd Umklapp 
light cones. X±zp at 0 near 60° has a transition dipole 
strength of ^ 0.0b2D at the main light cones, ^ O.OOII) 
at 1st Umklapp, and ^ 0.00051) at 2nd Umklapp light 
cones. D is the dipole strength of the intralayer excitons 
[401150]. Eq. 0 shows that the strength and polarization 
of T>r'r at the six-fold light cones have no dependence on 
the layer translation tq which appears in a phase factor 

g-ZTK-ro 

Translation dependence at commensurate stacking .— 
Eor MoSe 2 /WSe 2 or M 0 S 2 /WS 2 heterobilayers, neglect¬ 
ing the ^ 0.1% lattice mismatch, they have commen¬ 
surate stacking at 0 = 0°, 60°, 21.8° or 38.2°. At 
0 = 0°,60°, all main and Umklapp light cones merge 
into a single cone at Q = 0. At 0 = 21.8°, 38.2°, the 
six 2nd Umklapp light cones merge into a single one at 
Q = 0 (c.f. Eig. [^d)). Eor the light cone at Q = 0, 
the optical transition dipole is the superposition of con¬ 
tributions from multiple (tk, t'k') pairs each of which is 
associated with a distinct phase factor (c.f. Eq. 

Q). Their interference then gives rise to the different fq 
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dependence for the aH- and cr— components of r m- 

er-T>r'r oc (5) 

e-r-T>r'r OC + e-*^(<?3K.ro+^) _,_ g-iT(C|K.ro+^) 

Thus both the strength and the elliptical polarization of 
the transition dipole change with fq [51], as shown in 
Fig. g The above relations dictate that at fq = 0 and 
^^ 3 ^^ the in-plane transition dipole of X^-'r becomes fully 
circularly polarized, while at fq = both e± com¬ 

ponents of l^r'T vanish, in agreement with the ab initio 
calculations in Fig. [54]. Only at these three fq, the 
commensurate heterobilayers have the 3-fold rotational 
symmetry. For other fq, this rotational symmetry is ab¬ 
sent and elliptically polarized emission by at Q = 0 
is allowed. For 0 = 0°, the transition dipole is found to 
be strongest ^ 0.16D at fq = 0 (i.e. A A stacking, c.f. 
Fig. [^a)), while for 0 = 60°, the transition dipole is 
strongest ^ 0.15D at fq = g stacking, c.f. 

Fig. gb)). 

Discussions .—For interlayer excitons relaxed to the 
bottom of the energy dispersion, their photon emission 
need the assistance of phonon or impurity scattering into 
the light cones at the finite velocities. The inefficiency 
of such phonon or impurity assisted radiative recombi¬ 
nation process, compared to the direct recombination of 
intralayer excitons, can be relevant to the long lifetime 
(exceeding nanosecond) observed in the time resolved 
photoluminescence (PL) [19], which can be interesting 
for future study. The observation at the same time im¬ 
plies the inefficiency of non-radiative recombination of 
interlayer excitons. The energy of the emitted photon is 
EiiCl) plus (minus) the phonon energy if phonon absorp¬ 
tion (emission) is involved. The distribution of exciton 
kinetic energy and the phonon energy can account for 
the large broadening of PL. Our results also explain why 
strong interlayer exciton PL is observed only at twist 
angle 0 near 0° or 60° m, where the main light cones 
get sufficiently close to the bottom of energy dispersion. 
There are other twist angles {0 near 21.8° or 38.2°) where 
the 2nd Umklapp light cones are close to Q = 0, how¬ 
ever, the weak dipole of these light cones renders the PL 
inefficient. 
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